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SUMMARY/OVERVIEW: 

The  computationally  accurate  and  efficient  prediction  of  turbulent  reactive  flows  remains  an 
important  research  topic  for  several  reasons,  one  of  which  is  that  the  modeling  of  chemistry  is 
still  cumbersome,  relatively  inaccurate,  and  cannot  portray  the  complex  reactions  associated  with 
real  fuels.  To  reduce  computational  costs  and  enable  the  description  of  concurrent  reactions  from 
complex  fuels  containing  a multitude  a species,  a new  concept  is  here  proposed  that  relies  on 
statistical  instead  of  deterministic  concepts  for  modeling  chemical  kinetics  of  oxidation 
reactions.  In  this  new  concept,  each  species  has  a time  scale  associated  with  it,  and  the  goal  is  to 
solve  only  for  the  significant  scales  of  the  problem  and  model  those  scales  that  are  not 
significant.  A set  of  reaction-species  mole  fractions  is  considered  a vector  and  the  objective  is  to 
map  this  vector  into  a much  smaller  multidimensional  vector  space  of  base  chemical  fragments. 
The  goal  is  to  inquire  whether  such  a base  can  be  found  for  oxidation  reactions,  and  whether  one 
can  derive  corresponding  reaction-coordinate  rates  for  the  base  that  portray  the  energetics  of  the 
system.  In  essence,  there  is  a functional  mapping  representing  the  transformation  from  the  mole 
fraction  PDF  of  the  species  to  the  mole  fraction  PDF  of  the  base  and  the  interest  is  not  in 
calculating  this  functional  but  its  moments  (e.g.  statistical  means,  variances,  etc.)  needed  to 
adequately  describe  the  fragment  vector. 

TECHNICAL  DISCUSSION 

The  first  step  in  this  new  research  was  to  concentrate  on  finding  a potential  set  of  base  fragments. 
The  transformation  between  the  vector  of  the  species  and  the  base  fragments  vector  would  be 
represented  by  a matrix  [wc]  having  elements  representing  the  decomposition  of  each  species  on 
the  base.  Then,  the  rate  of  change  of  the  base  species  is  related  to  the  rate  of  change  of  the 
reaction  coordinates  through  the  transformation  [wc][/J  where  the  matrix  [%]  has  for  elements  the 
stoichiometric  coefficients  of  the  initial  set  of  detailed  reactions.  An  example  of  a detailed  kinetic 
model  is  the  LLNL  set  of  reactions  for  n-heptane  oxidation  which  involves  a total  of  1 60  species 
and  1540  reactions.  The  requirement  for  this  reduced  set  of  reactions  for  the  base  fragments  is 
that  it  preserves  the  heat  release,  indicating  in  particular  that  the  heat  capacity  and  heats  of 
combustion  must  be  appropriately  modeled.  The  concept  is  that  although  there  is  a functional 
mapping  representing  the  transformation  from  the  mole  fraction  PDF  of  the  species  to  the  mole 
fraction  PDF  of  the  base,  if  one  wants  to  reduce  the  complexity  of  the  system,  it  is  not  necessary 
to  explicitly  know  this  functional;  the  goal  is  to  obtain  an  implicit  model  based  on  correlations 


61 


resulting  from  numerical  simulations  using  detailed  models,  because  this  implicit  model  would 
be  much  simpler  than  the  explicit  one.  The  following  describes  a proposed  base  that  is  relatively 
simple,  yet  potentially  accurate. 

The  guiding  idea  for  the  base  search  is  similar  to  that  of  group  additivity  [1];  however,  group 
additivity  accounts  not  only  for  interactions  of  a group  with  its  nearest  neighbor,  but  also  for 
interactions  with  groups  not  immediately  adjacent  and  for  steric  effects.  The  present  intention  is 
to  only  account  for  the  interaction  with  immediately  adjacent  groups  as  the  complexity  of  the 
correlations  in  [1]  must  be  decreased  if  computational  efficiency  must  be  improved.  Thus, 
although  the  idea  is  very  similar  to  the  ‘groups’  of  the  group  contribution  methods  of  estimating 
species  properties  [2],  the  present  focus  is  on  a combination  of  comprehensiveness,  accuracy  and 
relative  simplicity  in  that  only  first  order  (compositional)  effects  are  allowed.  Examination  of 
alkanes  and  olefins  with  a vinyl  radical  shows  that  the  linearity  of  thermal  properties  in  carbon 
(C)  number  implies  that  the  proper  building  blocks  are  ‘component  radicals’  listed  in  Table  1. 
Listed  in  Table  1 are  also  values  of  thermodynamic  properties  and  correlations  for  calculating  Cp, 
listed  here  as  NCP  to  indicate  that  these  are  components.  The  fitting  of  the  simple  formula  (much 
simpler  than  in  [1])  is  tested  by  computing  Cp/R  (R  is  the  universal  gas  constant)  and  is  plotted  in 
Figs. la  (alkanes)  and  lb  (alkenes);  the  agreement  with  the  correlation  is  remarkable  given  its 
simplicity. 

However,  for  species  with  double  bonds  not  at  the  end  of  the  chain  or  for  oxidized  species,  the 
simplicity  of  the  component  radicals  as  building  blocks  no  longer  rigorously,  but  only 
approximately,  holds.  To  calculate  combustion  energies  once  either  breakup  of  a hydrocarbon 
into  free  radicals  or  oxidation  reactions  occur,  the  assumption  is  made  that  one  can  associate  a 
constant  partial  energy  with  all  component  radicals,  independent  of  the  species  to  which  they 
belong.  The  data  [3-5]  show  that  although  this  assumption  is  not  entirely  correct  (when 
subtracting  contributions  from  CH2,  CH3,  etc.),  the  variation  across  species  is  not  that  substantial; 
an  average  over  species  is  taken.  The  error  in  the  overall  heat  of  formation  for  the  entire  molecule 
is  rather  small;  two-digit  accuracy  is  normally  reached.  In  Fig.  lc,  Cp/R  is  plotted  for  some  ethers 
using  the  information  of  Table  1;  although  there  are  minor  deviations  from  the  data,  mostly  at 
low  T,  the  fits  are  adequate.  The  lightest  species  are  not  subject  to  meaningful  decomposition. 
Free  radicals  occur  when  a species  splits  in  two  (e.g.,  hydrogen  abstraction)  by  a reaction, 
forming  species  with  unsatisfied  valence  bonds;  component  radicals  could  also  exist  as  free 
radicals.  Properties  of  light  species  and  free  (component)  radicals  are  given  in  Table  2.  The 
corrections  for  the  heavier  free  radicals  may  be  estimated  using  the  increments  in  hydrogen  bond 
strengths  as  given  in  [6].  Figure  Id  illustrates  Cp  fits  for  a few  free  radicals,  and  good 
reproduction  of  the  data  is  shown.  Similar  fits  (not  shown)  for  the  light  species  of  Table  2 were 
equally  successful. 

Thus,  it  appears  that  the  base  will  be  formed  of  1 3 component  radicals,  1 1 free  radicals  and  5 
light  species.  A total  of  29  fragments  will  compose  the  base,  having  excluded  the  final  products 
of  combustion  (i.e.  CO2  and  H2O),  O2  and  N2,  for  which  the  mole  fraction  can  be  computed  from 
atomic  conservation. 
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Radical 

8h° 

8hc 

8hf 

ah 

bh 

CH3  (methyl) 

-42.0 

714 

188 

3.137 

3.433 

CH2  (methylene) 

-20.8 

614.3 

411 

2.784 

2.812 

CH  (methylidyne) 

s-7 

2.068 

3.038 

C2H3  (vinyl) 

62.5 

1212 

237 

4.432 

4.468 

C2H2  (CH=CH)° 

73±2 

1102±2 

not  relevant 

3.82 

3.72 

C2  (C=C,  dicarbon) 

230 

1017 

608 

2.919 

1.796 

HC2  (HC=C,  ethynyl) 

227 

1135 

339  ** 

3.937 

2.2515 

CO  (keto) 

-133±2 

260.5 

23 

2.892 

1.531 

HCO  (formyl) 

-124+3 

390.3 

167 

3.37 

2.92 

HO  (hydroxyl) 

-166±6*** 

-45 

204 

2.084 

1.0945 

-205+1 16/n 

-84+1 16/n 

215-1 16/n 

4.194 

1.908 

OO  (peroxy)* 

-16-13n 

-16-13n 

not  relevant 

3.65 

2.135 

EUS  £ S SSMIHHHI 

-74-13n 

-74-13n 

323+13n 

-1.46 

2.44 

0 trans:  -2,  cis:  +2;  for  a molecule  with  one  or  more  vinyl  also  present  (i.e.,  polyene),  8Cp/R  is  increased 
by  12%;  **  Alternate  value  (op.  cit.  NIST)  is  250  from  older  data;  ***  value  = -180  if  combined  with  CH 
radical;  * ‘n’  is  the  number  of  carbon  atoms  in  attached  radicals;  ® effective  n = 1.16  with  ring  or 
branched  attached  radical;  * max  value  An  = 2 for  an  attached  chain,  effective  n = 3 for  a ring  with  n > 4. 
If  at  least  one  vinyl  also  present,  then  add  An  = +1.  For  a ring  with  n = 2 or  3,  a correction  5h°  « 90 
applies.  (5Cp/R  fit  is  for  ring  ethers;  for  a chain,  8Cp/R  = 2.5  );  Note:  if  the  keto  (or  formyl)  radical  is 
attached  to  oxygen  (or  hydroxyl),  then  a correction  8h°  ~ -90  applies. 

Table  1.  Thermodynamic  properties  of  component  radicals.  (8h°,  8hc,  8hf  heats  of  formation, 
combustion,  and  component-to-free  transition,  respectively,  in  kJ/mol;  constants  ah  and  bh  for 
partial  molar  heat  capacity  in  the  form  8Cp/R  = ah  + bh  ln(T/T0);  T0  = 25  °C  = 298.15  K) 


Molecule/radical 

h° 

hc 

ah 

bh 

h2 

0.0 

241.5 

3.282 

0.400 

02 

0.0 

0.0 

3.476 

0.5663 

n2 

0.0 

0.0 

3.388 

0.469 

c 

717 

1111 

2.50 

0.0 

h2o 

-241.5 

0.0 

3.688 

1.217 

co2 

-393.5 

0.0 

4.690 

1.390 

N 

473 

473 

2.50 

0.0 

H 

218.0 

339 

2.50 

0.0 

HO 

38 

159 

3.385 

0.3637 

HOO 

10.5 

131 

4.150 

1.307 

O 

249.2 

249.2 

2.536 

0.0 

CO 

-110.5 

283 

3.426 

0.4749 

HCO 

43.1 

558 

4.154 

1.2875 

CH4 

-74.6 

802 

3.797 

4.305 

ch3 

146 

902 

4.440 

2.249 

ch2 

390 

1025 

3.973 

1.3015 

63 


CH 

596 

1110 

3.220 

0.7136 

c2h3 

300 

1449 

5.1 

3.5 

hc2 

566 

1474 

4.434 

1.404 

c2 

838*  1 

1625 

4.58 

0.0 

NO 

90.3 

90.3 

3.533 

0.4508 

no2 

33.2 

33.2 

4.691 

1.151 

h2o2 

-136 

5.269 

1.880 

c2h2 

228 

1257 

5.368 

2.294 

C2H4 

52.5 

1323 

5.383 

4.676 

Alternate  value  of  832  from  CRC  tables 

Table  2.  Thermodynamic  properties  of  molecules  and  free  radicals.  (h°,  hc  heats  of  formation  and 
combustion,  respectively,  in  kJ/mol;  constants  ah  and  bh  for  molar  heat  capacity  in  the  form  C/R  = 
ah  + bh  ln(T/T0)  .) 


T (K)  b T (K) 


Figure  1.  Plot  of  Cp/R  vs  log(T)  where  T is  in  degrees  K.  Data  values  are  as  symbols;  lines  portray  the  fit 
constants,  a.  alkanes,  constants  from  Table  1;  b.  alkenes,  constants  from  Table  1;  c.  ring  ethers  (CH2)nO, 
constants  from  Table  1;  d.  some  free  radicals,  constants  from  Table  2. 
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